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The structural, electronic, mechanical and thermal properties of Iri-^Rh^ alloys were studied systematically 
using ab initio density functional theory at different concentrations (x = 0.00, 0.25, 0.50, 0.75, 1.00). A Special 
Quasirandom Structure method was used to make alloys having FCC structure with four atoms per unit cell. 
The ground state properties such as lattice constant and bulk modulus were calculated to find the equilibrium 
atomic position for stable alloys. The calculated ground state properties are in good agreement with the exper¬ 
imental and previously presented other theoretical data. The electronic band structure and density of states 
were calculated to study the electronic properties for these alloys at different concentrations. The electronic 
properties substantiate the metallic behavior of alloys. The first principle density functional perturbation theory 
as implemented in quasiharmonic approximation was used for the calculation of thermal properties. We have 
calculated the thermal properties such as Debye temperatures, vibration energy, entropy, constant-volume spe¬ 
cific heat and internal energy. The ab initio linear-response method was used to calculate phonon densities of 
states. 
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1. Introduction 

Platinum group metals are promising candidates for a wide range of applications. Iridium (Ir), Rho¬ 
dium (Rh) belong to this group and are the rarest elements found in the earth’s crust. Their incredible 
catalysis activity increases research interest in these rare elements lU. Iridium has high density, high 
melting temperature and the highest corrosion resistance. Its thermal stability at high temperatures and 
the highest shear modulus at room temperature made it a superior metal (SHU. It has been investigated 
for carbon materials protective coatings |5), heavy-metal-ion sensors (6) and Re-rocket thrusters m ap¬ 
plications. It is broadly used in organic chemistry for hydrogenation, hydroformulation, hydroboration, 
hydrosilylation, cy do additions (8]|9l and oxidation reactions for catalytic converters in automobiles. It 
remains a metal of high interest due to its unique properties for the scientific community (To][TTl during 
the past decades. Binary intermetallic alloys that contain a transition metal display interesting electronic, 
structural, optical and thermal properties. Recently, Ir- and Rh-base alloys are known as ultra-high tem¬ 
perature materials due to their high-melting temperatures, good high-temperature strengths and good 
oxidation resistances II211141 . Rh-base alloys have better oxidation resistance, lower density, lower ther¬ 
mal expansion coefficient and higher thermal conductivity than Ir-base alloys. These properties make 
Rh-base alloys more promising for ultra-high temperature gas turbine applications. To our knowledge, 
theoretical and experimental calculation are performed for Ir, Rh and their alloys with other metals but 
neither theoretical nor experimental calculations are available for Ir-Rh alloys. These alloys provide a 
class of systems exhibiting unique mechanical properties that make them attractive for structural ap¬ 
plications (15). In recent years, Ir-base super-alloys with good high-temperature properties have been 
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developed as promising candidates for ultrahigh-temperature applications to replace the traditional Ni- 
base super-alloys, which have limited temperature capabilities due to the rather low melting point of 
Ni ITgWTsI . These alloys are promising materials for high temperature and pressure applications, and 
currently they are being examined for use in diesel engine turbocharger rotors, high-temperature die 
and molds, hydroturbines, and cutting tools (191. Due to the improvements in simulation techniques in 
material science, it is now possible to study materials from bulk to nano scale devices (20). The use of 
first principles calculations offers one of the most powerful tools for carrying out theoretical studies of a 
number of important physical and chemical properties of condensed matter with great accuracy I21II22I . 
Electronic structure simulations based on density-functional theory (DFT) (23][24l have been instrumen¬ 
tal to this revolution, and their application has now spread outside a restricted core of researchers in 
condensed-matter theory and quantum chemistry, involving a vast community of end users with very 
diverse scientific backgrounds and research interests. The knowledge of thermal properties is essential 
for the study of mechanical properties (25). To the best of our knowledge, there are no experimental and 
other theoretical data for comparison, so we consider the present results as a predictive study for the first 
time, which still awaits an experimental confirmation. 


2. Computational method 

The first principles investigation was performed using a pseudo-potential plane wave (PP-PW) method 
as implemented in QUANTUMES ESPRESSO (26). Special Quasirandom Structure method proposed by 
Zunger and coworkers (32) was used to make FCC structure with four atoms per unit cell alloys. SQSs 
are specially designed for a small supercell that is computationally feasible for DFT calculations. Vander¬ 
bilt ultra soft pseudo-potential parameterized by Perdew and Zunger (27][28) was used for calculations. 
We used ultrasoft Vanderbilt formalism with local density approximation (EDA) (29) for the exchange 
correlation energy of electrons. The special /c-points were integrated to the sampled Brillouin zone us¬ 
ing the Monkhorst-Pack method (31). A special /c-points mesh of 14 x 14 x 14 was used to produce an 
irreducible Brillouin zone. Pseudo-wave functions were expanded in a plane wave basis set using the 
cut-off energy of 25 Ry for all concentrations. The chosen cut-off energy and /c-point mesh ensure conver¬ 
gence with an accuracy of 10~® Ry. Density Functional Perturbation Theory (DEPT) and Quasi Harmonic 
Approximation (QHA) code developed by Baroni et al. (26) were used to investigate the phonon density 
of states and thermal properties from optimized structure. The dynamical matrices at arbitrary wave 
vectors were obtained using the Fourier transformation-based interpolations to calculate thermal prop¬ 
erties. A12 X12 X12 ^7 point mesh was implemented to obtain dynamical matrices of force constants in the 
irreducible Brillouin zone at E point. The linear response approach method was used to obtain the curves 
of phonon density of states within the framework of DEPT (26) as implemented in Quantum ESPRESSO. 
The parameters such as cell dimension, energy cut-offs and the /c-points used in the study were obtained 
from their respective convergence tests. These convergence parameters provide a fast way toward an 
optimized structure. 


3. Results and discussion 

3.1. Structural properties 

We have found the ground state properties of Iri-^Rh^ alloys at x = 0, 0.25, 0.5, 0.75, and 1.00. The 
total energy per unit cell was computed at various lattice parameters to find the equilibrium lattice con¬ 
stant. The equilibrium lattice constant was found by minimizing the total energy of system. Murnaghan’s 
equation of state was used to evaluate the optimized lattice constants and bulk moduli (33). 

The structural parameters i.e., lattice constants and bulk modulus for pure Ir, Rh and their alloys were 
calculated and reported in tablealong with experimental and other theoretical results. The optimized 
lattice curves for Iri-^Rh^ alloys at x = 0, 0.25, 0.5, 0.75, and 1.00 are shown in figure[^ Figure[^shows 
the variation of lattice constants and bulk modulus for Iri-^Rh^ evaluated using Vegard’s law. It is seen 
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Table 1. Calculated lattice parameters and bulk modulus of Iri-^Rh^ alloys compared with experimental 
and other theoretical results at x = 0, 0.25, 0.5, 0.75, and 1.00. 


Composition 

Fattice constant (Angstrom) 

Bulk modulus (Gpa) 

X 

This work 

Experimental 

Other 

This work 

Experimental 

Other 




calculation 



calculation 

0.00 

3.831 

3.84“ 

3.88^ 

360 

355^ 

386? 

0.25 

3.81394 



337.375 



0.50 

3.80373 



315 



0.75 

3.79589 



294.5 



1.00 

3.792 

3.803‘' 

3.847'^ 

269.5 

268.7-^ 

259.6^ 


^ reference (M), ^ reference ^ reference (36), ^ reference (37). 

® reference (38), ^ reference IMl . ^ reference |40l, ^ reference l4ll . 




Concentration X Concentration X 

Figure 1. Calculated lattice parameters and bulk modulus in comparison with Vegard’s law at experimen¬ 
tal values for Iri-^Rh^ alloys: (a) lattice constant, (b) bulk modulus. 


that with the change of Rh concentration in Iri-^Rh^ alloys, the lattice constant and bulk modulus show 
a considerable change from Vegard’s law which may be due to the lattice mismatch in these alloys. 

The calculated results are compared with experimental and other theoretical results for pure Ir and 
Rh. However, the experimental or theoretical results for their alloys are not yet available in the literature 
for comparison. These calculations may be a prediction for researchers for future. The values of the cal¬ 
culated lattice constant and the bulk modulus for Iri-^Rh^ alloys turned out to decrease with an increase 
of the Rh concentration. 

3.2. Electronic properties 

Electronic properties depend on the electronic configuration in materials, particularly on the exis¬ 
tence of prohibited regions of energy and on the magnitude in their electronic excitation spectra. De¬ 
tailed band-structure calculations are needed to understand the electronic properties of any material. 
Thus, the band structures for Iri-^Rh^ along various high symmetry directions were calculated at equi¬ 
librium lattice constants. The electronic band structures are shown in figure The Fermi level for the 
band structure was set to be 0 eV. It can be seen from figure (a)-(e) that electronic bands for all con¬ 
centrations overlap at the Fermi level. The electronic structures calculated for different concentrations 
clearly show the metallic nature of the materials since the energy bands intersect at the Fermi level. 
Conductivity in Ir metal is dominated by overlapping of 5, p and d bands. 

For Rh metal, the conductivity is dominated by a narrow d-band with some empty states overlapped 
by a broad free electron 5-band. Their alloys show a metallic nature due to s, p, d of Ir atoms and s, d of 
Rh atoms bands overlapping. So, these alloys lead to archetypical transition metals, in which a narrow d- 
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Figure 2. Structural optimization plots for Iri_xRhx alloys: (a) at x = 0.00, (b) at x = 0.25, (c) at x = 0.50, 
(d) at X = 0.75 and (e) at x = 1.00 concentrations. 


band containing some empty states was overlapped by a broad free-electron like 5-band which dominated 
the conductivity 1^ . 

The calculated electronic band structures along the principal symmetry directions are presented in 
figure [^(a)-(e). The overall band structure shows a similar metallic nature. The main hybridization re¬ 
mains between d bands of both Ir and Rh when mixed in order to form their alloys. The Fermi energy 
decreases with an increase of Rh percentage in Iri-^Rh^ alloys. The calculated electronic band structure 
is in good agreement with the previous result for pure Ir and Rh (43]|44l. From the calculated electronic 
band structure for Iri-^Rh^ alloys, it can be seen that the band overlapping increases with Rh concentra¬ 
tion which causes a decrease of resistivity with Rh concentration |4^ . 

The density of states plays an essential role in studying many physical properties. Total energies of 
material can be calculated from the knowledge of the density of states. The number of electrons within 
the Fermi surface can be used to determine the nature of materials. The Fermi energy is determined by 
the Density of States (DOS) which provides information on the Fermi energy level. To further study the 
nature of electronic band structure, we have also calculated the total and partial DOS for these alloys for 
different concentrations at an ambient pressure and presented them in figure[^(a)-(e). In figure[^(a), for 
pure Ir, the lowest lying bands were due to s, p states and the higher energy states were mainly due to 
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d states, and the conductivity was dominated by 5, p and d hybridization states. In case of Iro. 75 Rho .25 
in figure [^(b), the bottom state of the valence band was dominated with 5 , p states of Ir and Rh 5 state. 
The upper valance band region was dominated with Ir d state with some contribution of other states. 
In figure[^(c), with Iro. 50 Rho. 50 , the lowest part of the valance band was mostly contributed by 5 state of 
both metals and Ir p state but the higher energy level was dominated by a broad contribution of d-states 
of these metals with little part of other states. In figure|^(d), for Iro. 25 Rho. 75 , 5 , p states of Ir and Rh 5 state 
for this alloys dominated the bottom of the valence band. In the upper band, Rh d state was much wider 
than the Ir d states. In figure |^(e), with pure Rh, the bottom valence levels were mostly dominated by 5 
states but the higher valence band was dominated with d states. 



Figure 3. Calculated band structures for Iri-^Rh^ alloys: (a) at x = 0.00, (b) at x = 0.25, (c) at x = 0.50, (d) 
at X = 0.75 and (e) at x = 1.00 concentrations. 


From the analysis of the band structure of Iri-^Rh^ alloy, it was observed that, at different concen¬ 
trations of these metals, the band structure shows a similar behavior that differs from each other mainly 
by the energy level of each band relative to the Fermi level. In pure Ir, the conductivity was mostly 
dominated by s, p and d hybridization but for pure Rh, the conductivity was prominent with s and d 
overlapping. The electronic energy density of states are compared with the previously reported result 
for cubic Ir and Rh (46]|47l. When these metals were mixed to make alloys, strong hybridization between 
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Figure 4. Density of states for Iri-^Rh^ alloys: (a) at x = 0.00, (b) at x = 0.25, (c) at x = 0.50, (d) at x = 0.75, 
and (e) at x = 1.00 concentrations. 


the Ir and Rh d states with small contribution of s states dominated the conductivity. Ir has a lower elec¬ 
tronic DOS than Rh at the Fermi level and it increases with an increase of Rh concentration in the alloy. 
The material having a lower electronic DOS at the Fermi level often characterizes a more stable structure 
(48]Us). Thus, Ir is a more stable structure than Rh and it decreases with Rh concentration. 
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3.3. Thermal properties 

Thermal properties of metals are of interest below the melting temperature. We have used a quasi¬ 
harmonic approximation and phonon density of states to obtain thermal properties. Helmholtz free en¬ 
ergy AF, internal energy AE, entropy S and constant-volume specific heat Cy, at zero pressure are calcu¬ 
lated. Phonon DOS has a strong impact on thermodynamic properties. The electron excitation becomes 





0 50 100 150 200 250 

w/ cm'^ 


Figure 5. Phonon density of states for Iri_xRhx alloys: (a) at x = 0.00, (b) at x = 0.25, (c) at x = 0.50, (d) at 
X = 0.75, and (e) at x = 1.00 concentrations. 

easier with higher phonon DOS (50). Figure [^depicts the phonon density of states of Iri-^Rh^ alloys with 
different concentrations. The total phonon DOS for pure Ir was shown in figure[^(a) and has a maximum 
DOS at frequency 112.5 cm~^. In figure |^(b), two higher peaks in phonon DOS occur, first peak is dom¬ 
inated by Rh at frequency of 115.5 cm“”and the second was due to Ir atoms at 162 cm“^ frequency. In 
figurejsjfc), three high peaks for phonon density of states occurs at a frequency of 103.5 cm“\ 126 cm“\ 
164.25 cm“^ which are due to movements of Rh, Ir and Rh and Ir atoms, respectively. In figure|^(d), the 
contribution of total phonon DOS was dominated by Ir and Rh atoms at low frequency of 108 cm“^ Con- 
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versely, phonon density of states is mainly dominated by Rh atoms after this frequency. In figure (e), 
for pure Rh, two higher peaks occur in phonon DOS, the first peak is dominated mainly at 155.25 cm“^ 
frequency and the second peak with frequency 166.5 cm“^. 



Figure 6. Debye temperature variation with tem¬ 
perature for Iri-xRhx alloys at x = 0.00,0.25,0.50, 
0.75,1.00. 



Figure 7. Vibration energy variation with temper¬ 
ature for Iri_xRhx alloys at x = 0.00, 0.25, 0.50, 
0.75, and 1.00. 




Figure 8. Entropy variation with temperature for Figure 9. Specific heat at a constant volume vari- 

Iri_xRhxalloys at x = 0.00,0.25,0.5,0.75, and 1.00. ation with temperature for Iri-^Rhx alloys at x = 

0.00, 0.25, 0.50, 0.75, and 1.00. 

To further study the thermal properties, the Debye temperatures, vibration energy, entropy, constant- 
volume specific heat and internal energy are plotted in figures [6}fTo| for this alloy. Debye temperature 
determines the thermal characteristics of materials and it is closely related to many physical properties 
such as specific heat (51). The materials with high Debye temperature are associated with higher thermal 
conductivity. The knowledge of thermal conductivity and melting temperature is essential for develop¬ 
ing and manufacturing electronic devices (52|[53). The Debye temperature plays an important role in 
the field of thermo-physical properties of materials. Debye temperature was maximum at 0-20 K, then 
it decreases rapidly and approaches a constant for higher temperature values. Entropy, constant-volume 
specific heat and internal energy values increase quickly at a lower temperature and becomes converged 
to constant values with high temperature as in figureVibration energy remains constant for low tem¬ 
perature and decreases inversely as temperature increases. Entropy, constant-volume specific heat de- 
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creases with Rh concentration but increases for Iro. 50 Rho. 50 , then continues the decreasing behavior for 
further Rh concentration. The increasing behavior may be due to the maximum lattice mismatch at 50 
percent concentration and these quantities decrease with the further increase in Rh amount. It is seen 
that the internal energy and vibration energy increases slightly with Rh concentration but decreases for 
Iro. 5 oRho. 5 o concentration and then continues to progress. From figurej^ it can also be seen that Rh has 
a greater vibration energy compared to Ir. However, entropy, constant-volume specific heat, internal 
energy, vibration energy and Debye temperature show a different behavior at 50 percent concentration 
which may be due to a greater lattice mismatch. High entropy possesses many attractive properties, such 
as high hardness (54), outstanding wear resistance (55), good fatigue resistance characteristics (56), ex¬ 
cellent high-temperature strength (53, good thermal stability (58) and, in general, good oxidation and 
corrosion resistance (59). These properties suggest a great potential in a wide variety of applications. 



Figure 10. Internal energy variation with temperature for Iri-^Rh^ alloys at x = 0.00, 0.25, 0.50, 0.75, and 
1 . 00 . 

Our calculated results of thermal properties for pure Ir and Rh are in good agreement with the pre¬ 
viously published data (60lf63 . There were no experimental or theoretical results for their alloys. There¬ 
fore, the results reported in this work will be predictive for future. From the above investigations of 
thermodynamics properties, it can be demonstrated that Iri-^Rh^ with x = 0.5 has a larger entropy and 
constant-volume specific heat due to a maximum lattice mismatch within the considered range of tem¬ 
peratures but the vibration energy and the Debye temperature were minimum for x = 0.5. 


4 . Conclusions 

The first principle method was employed to analyze the structural, electronic, mechanical, and ther¬ 
mal properties of the Iri-^Rh^ alloys with a four atom unit cell. The pseudopotential scheme was used 
to study the alloys at different concentrations (x = 0.00, 0.25, 0.50, 0.75,1.00). The following conclusions 
were drawn: 

1. For the Iri-^Rh^ alloys, by increasing Rh concentration, the lattice constant and the bulk modu¬ 
lus decrease. The variations in the calculated lattice constant and the bulk modulus are slightly 
different from Vegard’s law. 

2. In electronic properties, the electronic bands overlap at the Fermi level, and the Fermi energy 
decreases with an increase in Rh concentration but the band overlap increases. This is because 
Rh has a larger electronic conductivity than Ir. In pure Ir, a narrow d state was overlapped by a 
broad free-electron s and p but pure Rh has an overlap between 5 and d only which dominates the 
conductivity in these metals. When Ir and Rh were mixed to form alloys, the main hybridization 
was between d states with some contribution from other states. 
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3. Thermodynamic properties such as phonon density of states, the Helmholtz free energy, the phonon 
contribution to the internal energy, the constant-volume specific heat, and entropy were calculated 
using the quasi-harmonic approximation. The minimum value of phonon density of states is for 
the alloy having Ir 0 . 50 Rh 0.50 concentration. Debye temperature was maximum at 0-20 K, then it 
decreases rapidly and approaches a constant for higher temperature values. Entropy, constant- 
volume specific heat and internal energy values increase quickly at a lower temperature and be¬ 
come converged to constant values with high temperature. Entropy, constant-volume specific heat 
decreases with Rh concentration but increases for Iro. 50 Rho. 50 , then continues a decreasing behav¬ 
ior for further Rh concentration. Vibration energy remains constant for low temperature and de¬ 
creases inversely as temperature increases. It is seen that the internal energy and the vibration 
energy increase slightly with Rh concentration but decrease for Iro. 50 Rho .50 concentration and then 
continue to progress. 
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Ab initio AocyiiA>KeHH5i CTpyKTypni/ix e/ieKTpoHHi/ix i Ten/ioBMX 
ByiacTi/iBOCTCM cn/iaBiB alloys 

HI. AxMefl, M. 3a(|)ap, M. LUaKi/i, M.A. Moy/pKapi 

/la6opaTopm KOMn'iorepHoro MOAe;iK)BaHH?i, (t)i3i/iMHi/iCi (t)aKy/ibTeT, Ic;iaMCbKi/iCi yHieepcuTeT m. Baxaea/inyp, 

M. Baxaea/inyp 63100, flaKi/icTaH 

CrpyKiypHi, e/ieKTpoHHi, MexaHiMHi i ren/ioBi B/iacTHBOcri cn/iaBiB Iri-^Rh^ ci/iCTeMari/iMHO AOC/iiAxeHO 3 bi/i- 
KopucTaHHfiM ab initio leopiT (t)yHKi4ioHa/iAy rycTMHi/i npi/i pi3Hi/ix KOHi^eHTpai^ifix (x = 0.00, 0.25, 0.50, 0.75, 
1.00). Byyio BUKopucraHO MeroA cnei^ia/ibHoi KBaaiBi/inaAKOBOi CTpyKiypn MOAe/iK)BaHHfi cn/iaBiB, u^o Ma- 
KDTb CTpyKTypy FCC 3 MOTupivia aroMaMi/i Ha OAHHi/iMHy KOMipxy. By/ii/i oOni/ic/ieHi B/iacTHBOcri ocHOBHoro craHy, 
Taxi 9\K nocTiwHa r'paTKi/i ra o6'eMHi/iM Mo^y/ib npyxHocri, /\n9] Toro, 3HaHTi/i piBHOBaxHi no/ioxeHH?i aroMiB 
fl,n9\ cTiHKHx cn/iaBiB. 06Mnc;ieHi ByiacTHBocii ocHOBHoro craHy Ao6pe y3roA>KyK)TbCfi 3 eKcnepi/iMeHTa;ibHi/iMi/i 
ra iHLUHMH paniLue orpHMaHi/iMH TeopeTHMHHMH A^hi/imh. 3 MeroK) BHBMenn^i e/ieKrpoHHi/ix B/iacTHBocreCi 141/ix 
cn/iaBiB npH pianiix KOHi4eHTpai4mx o6Mi/ic;ieHO e/iexipoHny aonny crpyKrypy ra rycTHHy craniB. E/ieKrpoHHi 
B/iacTi/iBocri oOrpyHTOByioTb Mera/iiMHy noBeAiHxy cnyiaBiB. J\i\^ o6Mi/ic;ieHH5i ren/ioBHX xapaKrepi/iCTHK 6y;ia 
Bi/iKopi/icrana nepiuonpi/iHi^i/inHa reopm sOypeHb (t)yHKi4ioHa/iy rycri/iHi/i, iMn/ieMeHTOBana y KBaairapMoniMHe 
Ha6/ii/ixeHHfi. HaMH o6Mi/ic;ieHo raxi ren/ioBi xapaxTepi/icri/ixH, fix reMneparypa /],e6afi, eneprm xo;ii/iBaHb, en- 
rponm, nHTOMa renyiocMHicrb npi/i nocriCiHOMy oO'cMi ra BHyrpiLUHfi eneprifi. By/10 Bi/ixopi/icrano ab initio ivieroA 
/liHiHHoro BiAryxyA/i^ o6Mi/ic/ieHHfi rycri/iHH craniB (t)OHOHiB. 

K/iiOHOBi c/ioea: e/ieKTpoHHi, craKTypHi i rennosi B/iacmBocri Mera/iiB n/iamHOBo'f rpyni/i 
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